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GENERAL  OBJECTIVES  0?  TASK  R 


Kar^y  of  the  vdde  variety  of  problems  associated  with  tlie  use  of  materials 
at  high  temperatures  occur  in  connection  with  advanced  propulsion  systems,  and  of 
these,  some  of  the  most  critical  and  complex  are  encountered  in  high  performance 
rocket  motors.  Present  trends  in  rocket  design  and  propellant  formulation  re¬ 
sulting  in  gas  flows  of  increased  temperature,  pressure,  and  corrosiveness  may 
be  expected  to  aggravate  the  materials  situation.  It  seems  self-evident  that 
future,  long-range  solutions  to  these  problems  must  rely  on  more  sophisticated 
approaches  and  broader  knowledge  of  the  behavior  of  materials  in  high  temperature 
gas  environments  than  is  characteristic  of  the  usual  "quick-fix",  or  expensive 
cut-and-try  test  procedures.  It  also  seems  likely  that  improvement  in  structural 
materials  themselves  (i.e.,  higher  strength  at  high  temperature,  higher  melting 
point,  etc.)  has  reached  the  point  of  diminishing  returns,  so  that  various  other 
subterfuges  must  be  tried. 

Research  performed  under  Task  R  al  the  Applied  Physics  Laboratory  or  its 
subcontracting  agfncies  is  intended  to  provide  some  of  the  fundamental  knowledge 
necessary  for  a  rational  understanding  of  the  behavior  of  materials  at  high 
temperatures.  This  is  a  very  broad  and  very  complex  field  involving  many  dif¬ 
ferent  scientific  disciplines,  V/hile  no  attempt  is  made  to  rigidly  limit  the 
scope  of  Task  R,  the  general  emphasis  is  on  appropriate  research  in  the  flow 
and  physical  chemistry  of  high  temperature  gases  such  as  are  characteristic  of 
advanced  solid  propellant  rocket  motors,  and  the  phenomena  basic  to  heat  transfer 
and  cooling  techniques  in  such  environments,  rather  than  in  the  properties  of 
materials  themselves. 


A,  A,  Westenberg 
Program  Coordinator 


SWIKARY 


I,  High  Tempera turo  Chemical  Kinetics  in  Laminar  Flames  . 

Tile  approach  taken  in  this  project  is  to  obtain  high 
temperature  chemical  kinetic  information  by  means  of  the 
detailed  analysis  of  laminar  flame  structure.  Net  rates  of 
appearance  or  disappearance  for  all  the  stable  species  found 
in  the  lAo  atmosphere  CH|,-'02  flame  are  presented  as  profiles 
through  the  flame  zone.  The  expected  relations  between  these 
rates  are  given  and  shown  to  be  borne  out  by  the  data.  The 
over-all  rate  constant  for  the  oxidation  of  CO  is  evaluated 
in  terms  of  an  empirical  mechanism  and  shown  to  be  essentially 
temperature  -  independent  in  the  range  1650-1950°K,  Comparisons 
with  0 tiler  work  are  made, 

II,  Thermal  Conductivity  of  Gases  by  a  New  Technique  . 

Extension  of  a  method  previously  developed  at  APL  for 
the  accurate  measurement  of  high  temperature  molecular 
diffusion  coefficients  to  the  measurement  of  thermal 
conductivity  is  being  attempted.  Improvements  in  both 
the  line  heat  source  and  the  wake  temperature  measuring 
equipment  are  described.  Previously  noted  difficulties 
with  the  perturbation  caused  by  the  velocity  wake  of  the 
source  wire  appe^ir  to  have  been  solved  satisfactorily. 

These  wake  corrections  are  described  and  results  on  room 
temperature  N2  and  He  of  considerable  precision  are  given. 

These  compare  well  with  values  measured  by  the  conventional 
technique. 

III,  Transport  Property  Studies  in  Dissociated  Gases 

Methods  of  measuring  diffusion  coefficients  and  thermal 
conductivity  in  gases  composed  of  labile  atoms  and  their  present 
molecules  are  being  explored,  A  proposed  diffusion  experiment 
is  described,  as  are  some  preliminary  experiments  being  per¬ 
formed  to  gain  experience  in  generating  and  detecting 
dissociated  atoms, 

IV,  High  Temperature  Chemical  Kinetics  from  a  Reacting  Point  Source  • 

The  idea  in  this  work  is  to  derive  chemical  kinetic 
data  by  analyzing  the  concentration  wake  downstream  of  a 
point  source  trace  gas  injected  into  a  laminar  flow  of  a 
second  gas  with  which  it  reacts.  The  theory  of  the  technique 
is  discussed,  and  the  apparatus  being  tried  out  is  described. 
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V,  Rocket  Kozale  Fluid  Dynamics  . 

This  project  aims  at  careAil  experimental  measurements 
of  the  flow  properties  in  a  rocket  nozzle  employing  a  typical 
solid  nropellant,  Pitot  pressure  surveys  in  nozzles  of  12,5® 
and  26^  total  divergence  angle  of  stations  corresponding  to 
tlic  same  expansion  ratio  show  appreciable  differences  which 
are  attributed  to  chemical  non-equilibrium  effects.  Some 
preliminary  results  with  a  fasvscanning  infrared  spectrometer 
to  measure  relative  gas  composition  in  the  nozzles  are  given. 
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I.  HIG{i  TEHFERATURE  CHSi^gCAL  Kir^CS  IK  lAHINAR  FLAKES 
(R.  H.  Fristronj  C.  Grunfelder,  S.  Favin,  A.  A.  Vfestenberg) 

Objective 

The  effect  of  chemical  reactions  is  one  of  tiie  important  problems  in  the 
dynamics  of  high  temperature,  high  speed  gas  flow.  Chemistry  and  aorod.v'nanics 
(and  sometimes  molecular  transport  processes)  are  strongly  coupled  in  various 
nozzle  flow  and  boundary  layer  i^enomena.  Advances  in  understanding  in  these 
areas  cannot  be  made  without  a  great  deal  more  knowledge  of  high  temperature 
chemical  kinetics  than  is  presently  available.  The  inroortant  reactions  are 
extremely  rapid  and  lie  in  temperature  regions  where  materials  present  problems. 
Therefore,  conventional  experimental  techniques  are  not  generally  applicable^ 
and  the  extrapolation  of  information  from  lower  temperature  reaction  rate  studies 
is  not  satisfactory.  Several  experimental  techniques  have  been  devised  to  sur¬ 
mount  these  difficulties.  Shock  tube  studies  are  useful  for  very  simple  reaction 
sj^tems.  Molecular  beam  studies  of  chemical  kinetics  (inelastic  scattering)  are 
in  their  infancy ,  but  offer  interesting  and  exciting  possibilities.  The  approach 
taken  in  this  project  is  to  study  high  temperature  kinetics  by  a  detailed  quanti¬ 
tative  analysis  of  laminar  flame  structure.  From  this  analysis  it  is  hoped  that 
eventually  the  kinetics  of  elemsntary  combustion  reactions  can  be  obtained.  The 
present  work  allows  the  determination  of  the  rates  of  appearance  or  disappearance 
of  all  of  the  stable  species  in  a  flame  as  a  function  of  position  (or  tirae)s 
tenperature,  and  chemical  composition.  This  is  the  same  type  of  data  obtained 
in  conventional  kinetic  studies,  and  a  much  wider  temperature  range  is  accessible. 


^  „ 


Exporipsntal  Techniques 

Tho  exporisental  techniques  used  to  obtain  the  teniperature,  con^sition^ 
and  aerod^nnanic  profiles  through  laninar  florae  zones  are  described  in  Quarterly 
Progress  Report  h'o,.  1  (Ref,  1). 


Kinetic  Analysis  of  tho  l/lO  Atmosphere  CH{^-02  Flarae 

Tho  set  of  data  obtained  in  a  flat  flame  (initial  mole  fractions i 

Xpu  «  0,0735 f  *  0,911j)  at  lAO  atmosphere  pressure  were  presented  in  earlier 

Quarterly  Progress  Reports  (Refs,  1^2)^  as  were  the  analyses  of  tliese  data  in 
terms  of  conservation  of  matter  and  energy.  During  the  present  report  period^ 
emphasis  has  been  on  interpreting  the  data  in  terms  of  chemical  kinetics. 


It  has  been  convenient  to  use  the  quantity  defined  by  the  relation 


G. 
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where  f^  is  the  mass  fraction  of  species  i  at  any  point  in  the  flame,  the 
diffusion  velocity,  and  v  the  macroscopic  gas  velocity  (mass  average  velocity). 

G^  thus  represents  the  fraction  of  the  total  mass  flux  at  any  point  due  to  species 
i,  and  includes  the  effect  of  diffusion.  It  may  then  be  shown  from  elementary 
species  continuity  considerations  that  the  voliunetric  net  rate  of  appearance 
or  disappearance  of  i  due  to  chemical  reaction  is  given  by 


£y\  ^  I 
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vriiere  p  is  gas  density.  is  molecular  weight,  and  z  the  distance  coordinate 
normal  to  the  flame  front.  When  the  usual  cgs  units  are  used,  has  the  units 
of  moles  cm  sec  . 
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The  rate  variable  was  calculated  for  all  stable  species  at  regular 
intervals  throughout  the  flane  by  roans  of  Eq.  (2),  uslr^  the  previously  derived 
profiles.  The  differentiation  was  done  nurorically.  Since  the  evaluation  of 
the  required  a  nuraerlcal  differentiation  of  the  original  conpositlon  profiles^ 
the  detemination  of  the  involves  Implicitly  a  second  derivative  of  the  experi- 
rontal  data.  This  is  an  obvious  and  unavoidable  source  of  error.  Since  errors 
my  be  expected  to  rou(^y  double  upon  each  differentiation^  this  whole  procedure 
probably  multiplies  original  errors  in  composition  by  a  factor  of  about  four. 

It  is  estimated  that  the  precision  (reproducibility)  of  the  composition  data 
is  about  2«3JSj  so  that  precision  in  the  after  the  data  analysis  may  be  about 
IC^,  This  says  nothing  about  other  sources  of  error,  which  are  discussed  in 
detail  elsewhere  (Ref.  3).  The  profiles  derived  in  this  way  are  given  in 
Fig.  I"l,  The  profile  for  Og  is  not  included  since  it  is  not  an  independently 
derived  quantity,  while  that  for  HCKO  is  not  shown  because  it  is  quantitatively 
negligible  compared  to  the  others.  The  rates  are  not  plotted  to  zero  value  at 
the  cold  boundary  side  of  the  figure  because  reliable  derivatives  of  the  G^  curves 
could  not  be  obtained  from  the  data  all  the  way  to  their  asyr?)totic  values,  Tlie 
rates  approach  zero  rapidly  to  the  left  of  the  first  plotted  points  of  the  curves, 
however.  It  should  be  noted  that  the  free  radicals  actually  present  in  the  flame 
are  not  analyzed  for  by  present  techniques,  and  show  up  in  the  stable  species 
concentrations.  Thus  their  influence  on  the  rate  profiles  cannot  be  assessed 
quantitatively  as  yet. 
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Socw  intorosUng  qualitative  observaiicrs  can  bo  rado  fron  those  rate 
profiles  first  of  aU,  It  will  be  not4Qd  Umt  the  flano  can  bo  convordontly 
divided  into  three  tones.  ‘Hm!  first>  extondii^  frcn  the  screen  (&  •  0)  to 
approxlnately  t  •  0,2  cn,  corresponds  to  a  region  wiiero  all  duinges  in  conpo- 
sition  and  tenperature  are  duo  sinply  to  transp<»rt  processes.  There  is 
essentially  no  chenical  reaction  until  a  tenperature  of  about  13<XPis  reached. 

The  second  tone,  covering  the  approxlnato  range  z  ■  0,2  -  0,3  cn,  is  where  nest 
of  the  fast  reactions  occur.  Here  the  CHj^  largely  disappears  with  fornation  of 
most  of  the  CO^  and  The  third  zone,  from  z  *0.3  cm  on  out  to  the  hot 
boundary  at  z^=^0.9  cm  (far  beyond  the  limit  shown  on  Fig,  I-l),  is  a  slow  re¬ 
action  zone  where  CO  is  converted  to  CO^.  This  is  the  CO  "afterburning”  region 
previously  noted  by  others  (Kef,  h). 

From  a  quantitative  viewpoint,,  the  data  may  be  used  to  derive  considerable 
information.  These  are  noted  as  follows: 

i)  V/ithout  specifying  any  details  of  the  actual  mechanism  whidi  may  pre¬ 
vail  in  the  flame  the  various  rates  should  be  related  to  one  another 
in  a  predictable  way.  Since  atoms  of  a  particular  element  are  neither 
created  nor  destroyed  in  the  combustion  process,  it  is  possible  to  write 
in  general  that 

K  -0  (3) 

i 

where  n^  is  the  number  of  atoms  of  an  element  in  a  molecule  of  species  i, 
Eq.  (3)  should  be  obeyed  for  each  element  throughout  the  flame.  Thus  for 

*^C0  "  “  ”  ^C02  -  ^CH20 


carbon 
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Figuro  I«-3  shows  a  corqiarison  of  tltoso  experirental  quantities^  i.e.; 

and  -  Kqq^),  The  contribution  of  was  neglected.  The 

agreement  is  generally  very  good.  TIu}  two  curves  have  similar  shapes, 
have  their  naxina  and  nlrdna  at  about  the  same  spatial  positions,  and 
the  absolute  nagnitudos  agree  quite  well  especially  in  the  long  CO 
aftorburnii^  region  where  has  vanished.  Similarly,  for  hydrogen. 


Eq.  (3)  yields 


*^H20  *  '  ■“  ^CHgO 


(5) 


and  the  comparison  of  the  two  sides  of  Eq,  (5)  is  given  in  Fig.  1-3 
(neglecting  q).  Again  the  agreement  is  generally  good.  The  worst 
deviations  occur,  as  usual,  in  the  region  of  most  rapid  reaction,  which 
right  be  expected.  The  fact  that  the  requirements  of  element  conservation 
are  met  reasonably  well  by  the  data  shows  chat  they  are  internally  con¬ 
sistent  and  the  analytical  treatment  used  to  derive  the  rates  is  sound, 
ii)  For  use  in  practical  calculations  (such  as  in  the  case  of  nonequilibrium 
nozzle  flow),  it  is  often  necessary  to  have  data  on  the  over-all  kinetics 
of  a  given  reactant.  Even  in  the  rare  case  where  a  reasonable  kinetic 
mechanism  is  available,  it  is  usually  so  complex  as  to  be  impractical 
to  use  in  calculations  which  are  already  very  involved.  From  the  over¬ 
all  (also  sometimes  denoted  by  the  execrable  term  “global”)  viewpoint, 
nothing  vrtiatever  is  implied  about  the  actual  elementary  reaction  steps 
in  the  process.  A  rate  expression  is  set  up  which  involves  the  concen¬ 
trations  of  the  reactants  and  an  adjustable  proportionality  constant 
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a  rate  constant  -  function  of  tereperature  in  general),  and 
tlie  experinental  data  are  fitted  to  this  expression.  It  is  strictly 
an  enpirical  procedure.,  justified  only  if  it  is  useful. 

The  ne thane -oxj’gon  flane  data  can  bo  treated  in  this  wai'.  The 
CO  oxidation  provides  a  goed  example.  Lot  us  write  the  expression 


dCCO^) 


r.  -  k  (co)(o,,)  -  k  ir  X, 


CO.. 


CO  0. 


(6) 


whore  the  parentheses  indicate  concentration  in  tnolos  cni*"^,  K  is  the 
total  molar  concentration,  and  the  X's  are  mole  fractions,  k  is  an 
empirical  second  order  rate  constant  liaving  units  of  cm^  raole**^  sec^^. 

We  presume  that  such  an  expression  holds  when  there  is  considerable  UgO 
present,  since  the  rate  of  oxidation  of  CO  with  very  little  H2O  is  strongly 
dependent  on  the  H^O  concentration.  It  is  better  to  write  Eq.  (5)  in 
terms  of  the  appearance  rate  of  CO2  rather  than  the  disappearance  rate 
of  CO5  since  the  latter  is  not  directly  determined  experimentally  except 
in  the  afterburning  region  (see  Eq.  (3))o  Using  the  experimental  values 
of  KpQ  ,  N,  and  the  mole  fractions,  the  rate  constant  k  was  determined 
throughout  the  region  where  there  were  valid  K..  data,  i.e.,  z  =  0.27 

UU2 

-  0,60  cm  corresponding  to  temperatures  l650-19li0°K,  The  values  of  k 
so  obtained  were  essentially  independent  of  temperature  and  centered 
around  the  value 


o  9  3  -1  ..1 

\i  -^0  X  10  cm  mole  sec 

The  temperature-independence  presumably  is  an  indication  that  the  CO 
oxidation  actually  is  a  free  radical  reaction  of  low  or  zero  activation 


energy 
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It  is  of  interest  to  conpare  Uus  value  with  those  obtained  by  other 
workers,  Feninore  and  Jones  (Hof.  5)^  using  flaroes  of  several  different 
hydrocarbons  at  atroosphoric  pressure,  concluded  that  the  rate  constant  defined 

9 

in  Eq.  (5)  had  a  value  of  about  1-2x10  in  a  sinilar  terperature  range.  They 
attributed  an  Arrhenius -t>T>e  tc*n}»raturc  dependence  to  k,  however,  and  assigned 
an  activation  enorg>’,  Friednan  and  C^Tphers  (Hof,  6),  using  a  propane-air  flame 
at  b6  mm  pressure,  obtained  data  wliich  can  be  corrected  to  conform  to  our 
definition  of  k.  The  result  at  i600°F»  was  about  LxlO^,  For  a  flame  of  CO 
and  O2  with  several  percent  of  H^O  present  burning  at  30  nan  pressure,  Friednan 
and  Nugent  (Ref.  7)  obtained  data  which  i'ields  a  value  of  about  IxlO^  at  i800°K, 
while  with  H2  in  the  original  mixture  burning  to  give  the  sane  amount  of  H^O 
the  values  of  k  computed  are  2-3  times  higher  than  this.  These  various  experi¬ 
mental  rate  constants  are  summarized  in  Table  I-l. 

Table  I  I.  Comparison  of  over-all  rate  constants  for  tlie  CO-O2  reaction  from 
work  of  various  authors.  Excess  H2O  present. 


Author 

Flame  Type 

k  (cm\ole"^sec 

This  work 

CHj^--02  at  76  nun 

8x10^  ( 1650-1 9h0°K) 

Fenimore  &  Jones 
(Refo  5) 

Hydrocarbon-air 
at  one  atm. 

1-2x10^  (1700-1900°K) 

Friedman  &  Cyphers 
(Refo  6) 

C^Hg  -air  at  ii6  mm 

IjxlO^  (l600°K)  (corrected) 

Friedman  &:  Nugent 
(Refo  7) 

CO  O2-H2O  at  30  mm 

1-3x10^  (1800°K) 
(corrected) 

All  of  these  data  were  obtained  by  flame  structural  studies.  No  other  data  by 
other  techniques  are  known  to  us,,  although  such  a  comparison  would  be  very 
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valuable.  In  any  case^  tlie  data  civen  in  Table  I-l  appear  in  reasonable  agree- 
raent  -  quite  ronarkable  agreement  coi^sidering  the  vide  variety  of  oxporir»ntal 
pressures.,  reactants,  and  tedmical  details  of  sa.'^pling.  analysis,  correction  for 
difftision,  etc.  Accordingly,  it  seece  likely  tJwt  the  over-cll  rate  constant 
for  the  reaction  defined  as  in  Eq,  f5)  *s  quite  reliably  established  in  Uie  indi¬ 
cated  temperature  range,  and  could  be  used  in  practical  computations  for  w!uch 
it  is  appropriate. 

Heat  Release  Rates  in  the  1/10  Atmosphere  0^  riame 

It  IS  of  some  interest  to  compute  the  volumetric  lieat  release  rate  profile 
from  the  experimental  data.  There  are  two  ways  of  doing  this,  which  can  then  be 
checked  against  each  other.  If  we  define  Q  as  the  rate  of  heat  evolution  (a 
positive  quantity).,  ther.  tne  simplest  way  to  compute  it  is  by  way  of  the  relation 

Q  —  I  \l  K,  (7) 

1 

where  is  the  absolute  molar  enthalpy  of  species  i.  The  negative  sign  is 
necessary  to  conform  to  the  above  sign  convention,  whlcl’.  is  contrary  to  the 
usual  thermodynamic  convention.  The  alternative  approach  is  to  use  the  energy 
conservation  and  species  continuity  equations  which  may  be  combined  to  give 


dT  V  Vi  Id  /..dTx 


(6) 


where  is  the  molar  heat  capacity  at  constant  pressure,  A  is  the  streamtube 
area  ratio  relative  to  some  reference  value,  and  is  the  gas  thermal  conduc¬ 
tivity. 
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Tho  rothod  of  conputing  Q  from  Eq,  (6)  involves  only  the  chemical  rates 
and  the  corresponding  enthalpies,  so  it  is  a  direct  neasure  of  the  rate  of 
enthalpy  change  per  unit  volume  due  to  reaction.  The  metliod  of  Eq.  (7)  regards 
Q  as  a  source  term  in  Uie  enthalpy  flux  conservation  equation  and  involves  the 
tempera ture  and  profiles.  TJui  results  of  both  approaclwjs  are  shown  plotted 
in  Pig.  I’-h.  Typically,  tho  worst  deviations  occur  in  tho  fast  reaction  zone. 
This  may  be  due  to  neglect  of  free  radical  effects,  since  presumably  if  those 
could  bo  sampled  and  their  rates  (positive  quantities  in  this  region)  de¬ 
termined,  they  would  contribute  positive  enthalpy  to  the  sura  in  Eq.  (6)  and 
thus  lower  the  value  of  Q  obtained,  bringing  it  more  nearly  in  agreement  with 
the  curve  determined  from  Eq.  (7).  The  two  curves  agree  very  well  in  the  CO 
afterburning  region. 

The  l/20  Atmosphere  Flame 

The  experimental  data  on  this  flame  have  now  been  obtained  and  are  being 
analyzea.  Results  oi  this  work  will  be  descril^d  in  the  next  Quarterly  Progress 
Report, 
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II.  TtgPJIAL  C0:n)lCTIVI7Y  Or  GASES  W  A  T5CHMIQUE 
(R,  E,  Walker^  }I.  doHaas,  A.  A.  Westenborg) 


Objective 

The  analysis  of  the  uroblen  of  coiwective  heat  transfer  both  in  the  lardnar 
and  turbulent  flow  rogiros  requires  tliat  the  properties  of  the  gases  involved  bo 
quantitatively  understood,  TI»o  coefficient  of  thernal  coixiuctivity  is  one  such 
property.  Heasurer/ent  of  thernal  conductivity  of  gases  has  not  been  a  forsaken 
field  of  rcsoardi.  and  several  coccson  gases  have  been  investigated  up  to  tenpera- 
tures  approaching  reaction  motor  temperatures.  The  primary  shortcomings  of  past 
and  present  research  in  the  study  of  thermal  conductivity  of  gases  for  application 
to  such  heat  transfer  problems  are?  a)  inability  to  apply  present  experimental 
methods  to  temperatiu'es  above  about  1(XX/^K  (most  measurements  are  made  at  sub¬ 
stantially  lower  temperatures  than  this),  b)  the  complexity  and  difficulty  in 
using  the  usual  apparatus  for  this  measurement  (thermal -conductivity  cells), 
and  c)  the  lack  of  measurements  on  species  (and  mixtures)  prevalent  in  reaction 
motors.  Other  shortcomings  could  be  listed. 

Recent  measurements  of  molecular  diffusion  coefficients  -  a  closely  re¬ 
lated  transport  property  (Ref,  1)  of  common  flame  gases  at  temperatures  in 
excess  of  1000°K  applied  a  simple  concept  (Ref,  2)  that  could,  in  principle, 
be  used  for  measuring  other  gas  transport  properties  (conductivity  and  viscosity). 
In  this  method,  a  point  source  of  a  trace  gas  is  located  in  the  center  of  a  uni¬ 
form,  heated  laminar  jet  of  a  second  gas,  A  series  of  gas  samples  removed 
immediately  dovnistream  of  the  source  and  quantitatively  analyzed  is  used  to 
determine  the  molecular  diffusion  coefficient.  The  technique  has  good  precision 
(+1-2/'3)  and  is  quite  adaptable  to  high  temperatures  and  various  gas  types. 
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Application  of  this  nethod  to  the  Reasurenant  of  Uiomal  conductivity  of 
gases  is  an  obvious  extension.  T}>e  source  of  trace  gas  is  replaced  with  a  point 
(or  lino)  source  of  licat  and  the  tenperatiire  rise  downstrean  of  the  heat  source 
is  neasiurod  with  an  appropriate  technique.  The  nathod  is  simple,  should  have 
good  precision  {♦2;4).  and  can  be  used  with  a  variety  of  gases  at  high  ambient 
temperature  (approaching  flame  temperatures).  Thus,  the  objective  of  this  re¬ 
search  is  to  employ  the  point-source  or  the  line*source  technique  in  the  maasurc- 
nent  of  gas  thermal  conductivity  and,  if  successful,  to  use  the  method  for 
obtaining  data  over  a  vdde  range  of  temperature  on  pure  gases  and  mixtures  of 
propulsive  interest. 

Theory  of  the  1‘lethod 

This  has  been  described  in  detail  in  Ref.  3. 

Experimental  Technique 

During  the  period  of  this  report  a  number  of  changes  in  the  experimental 
apparatus  have  been  made  from  that  originally  described  in  Quarterly  Progress 
Report  No.  1  (Ref.  h).  A  description  of  these  changes  follows, 

1.  Line  heat  source.  In  the  case  of  the  line  source,  heat  conduction  along 
the  wire  to  the  flow  tube  walls  where  it  is  attached  may  give  rise  to  a 
significant  error  in  the  determination  of  q,  the  heat  dissipation  per  unit 
length,  to  be  used  in  the  working  equation  (Ref.  3)  for  the  determination 
of  the  thermal  conductivity  of  the  gas.  The  proper  q  to  be  used  in  the 
working  equation  is  the  q  which  corresponds  to  the  central  portion  of  the 
source  wire.  The  temperature  decrease  at  the  wire  ends  due  to  heat  con¬ 
duction  to  the  v/alls  decreases  the  resistance  of  the  wire  in  this  region, 
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so  that  the  avoraco  q  for  the  entire  wire  differs  fron  the  q  in  the 
central  portion.  This  error  ra;/  be  reduced  by  using  a  wire  of  low 
thernal  conductivity  and  l(w  temperature  coefficient  of  resistivity, 
and  also  by  attaching  the  voltage  measuring  leads  to  the  saino  wire 
at  some  distance  from  the  walls.  In  view  of  these  facts,  a  line  source 
of  0,0005 -in, -diameter  platinum  -  10^  rhodium  wire  was  rx)unted  in  a 
3  cm  flow  tube  and  is  now  in  use.  This  source  wire  has  potential  leads 
of  tlie  same  wire  welded  to  it  at  points  approximately  h  m  from  the  walls, 

2,  Measurement  of  the  temperature  rise.  The  method  of  measurement  of  the 
temperature  rise  downstream  of  the  source  has  been  improved  by  the  use 
of  a  high  sensitivity  DC  breaker  amplifier.  The  output  of  the  dif¬ 
ferential  thermocouple  (Ref,  h)  is  applied  to  the  input  of  the  amplifier 
and  the  output  is  displayed  on  a  chart  recorder.  The  input  circuit  to 
the  amplifier  was  constructed  so  as  not  to  introduce  harmful  thermo¬ 
electric  potentials  and  was  shielded  to  reduce  the  effect  of  stray  fields, 
A  known  resistance  was  included  in  the  input  circuit  for  the  purpose  of 
applying  an  accurately  knovm  calibrating  voltage  to  the  amplifier.  Also, 
a  variable  negative  feedback  system  was  installed  in  the  amplifier. 
Adjustment  of  the  feedback  ratio  changes  the  effective  input  resistance, 
which  then  provides  a  check  on  Peltier  heating  of  the  thermocouples. 

The  advantages  of  the  above  system  over  that  of  using  a  potentiometer 
to  measure  the  differential  thermocouple  EMF  directly  are  brieflyj 

Beckman  Instruments,  Inc,,  Model  llj. 

A  Rubicon  High  Precision  Type  B  Potentiometer. 
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(a)  T:no  ni-oolraon  of  -jor.i lu-c  ^nt  in  'nrrc.nno^i  ;')\v,  \  \  lo 

apj*roximU»ly  ^1/2J  i.v.  This  nakes  ponsiblc  tlio  accui'ato  ncasuro- 
nsnt  of  snaller  tonroraturo  difforoncoa  and  tho  use  of  a  Uiorno- 
couplo  nysten  of  lover  (and  r.oro  precisely  known)  thcnwolectric 
power,  such  as  Pi-ra,  M:.  1 1  also  ailcvs  lower  soui'co  itoating 

rates  to  be  used. 

(b)  Tho  total  tine  required  to  rake  tho  nunb-jr  of  noasuromonts  nocos- 
sary  for  a  tliomal  conductivity  dotomination  at  a  given  temperature 
is  reduced  by  a  factor  of  3*  Tids  is  particularly  desirable  in 
regard  to  those  changes  in  the  entire  system  which  effect  the  indi¬ 
vidual  thermocouple  EI-F's  and  are  not  canceled  out  in  the  two 
opposing  thernocouples, 

(c)  Looking:  forward  to  measurements  at  higher  temperatures,  it  may  be 
expected  that  the  above-mentioned  uncompensated  drift  may  be  greatly 
enlarged  through  time  and  position  variations  in  the  gas  temperature. 

If  this  occurs,  it  may  become  necessary  to  chop  the  current  to  the 
wire,  thus  providing  a  square  wave  heat  source,  V/ith  the  thermocouples 
fixed  in  position,  adjustment  for  zero  vrould  be  made  during  the  time 

of  zero  current  and  a  single  potential  measurement  would  be  made 
while  the  current  is  flowing.  With  the  currently  alternately  off 
and  on  for  a  period  of  15  seconds  each,  it  is  then  required  that  the 
uncompensated  drift  be  small  for  approximately  this  length  of  time. 

The  present  amplifier-recorder  system  is  particularly  adapted  for 
this  method  of  measurement. 

The  reproducibility  in  the  thermal  wake  temperature  measurements  using 


the  line  source  and  instrumentation  described  above  is  better  than  +1/2% 
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Rcsultg  of  Hcasurera>nta  l‘girv:  tha  Lire  Soiurca 

In  Quarterly  Frogrosa  Report  I'o,  2  it  was  reported  that  Urn  raasured 
thomal  conductivity  of  nitrogen  at  usint:  tiio  line  sow'ce  varied  narkedly 

with  the  downs trearj  position  (x)  of  the  transverse  pi*ofile  rjoasurorents  (see 
Fig.  II-5  of  Ref,  3).  sluring  the  period  of  U»is  report  it  has  been  shown  that 
this  variation  of  neasurod  with  the  profile  x  is  due  to  non-unlfomity  in 
the  velocity  field  downstroan  of  the  wire  heat  source,  A  Uioorotical  study  of 
tlus  has  pointed  the  way  to  a  correction  to  the  data  which  has  given  neasured 
A's  Uu»t  are  independent  of  the  profile  x,  Tliis  will  be  shown  in  the  presentation 
of  results  for  nitrogen  and  hcliun  at  rcon  tenperalure,  obtained  by  using  the 
apparatus  described  in  the  previous  section. 

Table  II >1  lists  three  sets  of  ^M^asured  quantities  and  computed  thermal 
conductivities;  two  are  for  nitrogen  at  differing  rates  of  heat  release  and  the 
third  is  for  helium,  Sach  set  consists  of  raximum  profile  potentials, 
measured  at  a  numl>er  of  relative  downstream  distances,  x,  and  half  widths  of  the 
temperature  profiles,  z,  corresponding  to  a  (d-ven  E^/E  value  at  each  relative  x. 

In  the  section  labeled  "uncorrected"  there  aijpear  computed  A's  which  were  ob¬ 
tained  using  the  working  equation 


A 


2 


m 


(r-x) 


fE  r~' 

log 

-2  - 
E  '  r 

L  J 

where  a  (=  E/AT)  is  the  thermoelectric  power  of  the  downstream  thermocouple.  In 

2 

this  analysis  of  the  data  the  relative  x  were  displaced  to  make  xE^  independent 
of  X.  Fig,  II-l  includes  a  curve  of  the  uncorrected  A  plotted  against  the 
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1 

1  (to) 
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(cal/cnaQC®S 

Values 

^  X  10^ 

iu60 

26.60 

.955 

6.75 

>.oo 

23.65 

6.23 

6.60 

23.91 

1.123 

6.59 

3.00 

20.59 

6.3I2 

9.60 

19.91 

1.331 

6.1.1: 

11.00 

17.61 

6.31 

6.21;  (Rof.  5) 

13.60 

16.65 

1.513 

6.27 

15.00 

15.01; 

6.3I1 

17.60 

11:.  65 

1.717 

6.15 

19.00 

13.10 

6.25 

raan*  6,29 

8.35 

16.00 

1.263 

6.73 

iO.95 

13.07 

6.20 

11.35 

13.71: 

1.138 

6.57 

13.95 

11.59 

6,33 

15.35 

U.8I 

1.632 

6.L1 

17.95 

10.23 

6.32 

6.25  (Rof.5) 

20.35 

10.22 

1.838 

6.26 

22,95 

9.03 

6.30 

23.35 

9.51- 

1.91:7 

6.20 

25.95 

8.1.9 

6.28 

mean  -  6.3O 

3.35 

9.88 

1.169 

1:5.2 

5.85 

6.80 

36.9 

11.35 

8.78 

1.290 

1:3.1: 

6.85 

6.31: 

37.1 

6.35 

7.33 

1.190 

10.9 

8.85 

5.63 

36.9 

8.35 

6.32 

1.61:7 

39.8 

10.05 

5.06 

37.6 

35.9  (Ref.6) 

10.35 

5.69 

1.806 

38.1: 

12.85 

i:.67 

37.1 

112.35 

5.17 

1.92ii 

37.7 

U.85 

ii.33 

37.1 

ill. 35 

I1.76 

2.060 

37.3 

16.85 

1:.05 

37.3 

L— 

mean=  37.1 

^11  ^ 

A.  N^;  U  =  1]0  cm/sec,  q  =  9.12  x  lO"'’  cal  sec  cm""^,  T  =  300°K,  ~  *  2.00. 

B.  Np!  U  <=  IjO  cm/sec,  q  =  6.96  x  lO"^  cal  sec"^  cm"^,  T  =  301°K,  —  -  2.00. 


C.  He:  U  =  100  cm/sec,  q  =  8.75  x  lO”'^  cal  sec 


-1  -1 

cm 


T  =  298°K,  ~ 


f 


l.ij3. 
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profile  X  for  3«t  A.  "conroctecl**  section  of  Table  II-l  presents  results 
vJilch  were  obtained  bv  applying:  a  correction  vMch  was  necessitated  by  the  vel¬ 
ocity  wake  fron  tbe  heating  vlro.  The  treatnont  of  Uto  data  Including  the  wake 
correction  proceeds  as  follcvst 

(1)  Tiw  reasured  relative  x  are  displaced  to  nako  (r-x)/loc|~ 

Independent  of  x.  The  neu  downstrean  distances  are  designated  x^. 

(2)  TJie  reasured  k  are  adjusted  accordii^  to  E  ■  E  (1-A£  ),  Tite 

(o)2 

constant  A  is  obtained  by  requiring  x^E  to  be  equal  to  two 

c  R 

sera rated  x^  values, 
c 

(3)  r)  is  detemined  using 


fV 


2  2 
q  a 

*  '  • 


8n  jx  E 


(r-xc) 
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■»  E  V  r 


j  logi 

Table  II-l  lists  the  values  obtained  by  this  method  along  with  the  cor¬ 
responding  X  and  values.  It  is  seen  that  the  A  values  for  the  three  sets 
c  n 

of  data  are  no  longer  dependent  on  the  profile  x  value  and  that  the  average 
compares  favorably  \fith  the  measurements  of  others.  Fig,  II-l  shows  the  con¬ 
stancy  of  the  corrected  h  with  profile  x  for  data  sets  A  and  B, 

It  is  also  to  be  noted  that  the  origin  of  the  heat  source  is  nearly  the 
same  for  the  three  sets  of  data,  which  is  as  expected.  This  is  in  contrast  to 
the  uncorrected  data  vAiere  the  origins  differ  by  more  than  a  millimeter  among 
the  sets,  (Compare  the  x's  corresponding  to  =  11,00,  10,95,  and  10,85  of  the 
three  sets, ) 


In  almost  all  cases  it  is  sufficient  to  replace  this  quantity  by  (r-x)  and 
in  a  number  of  cases  z  is  small  enough  that  the  quantity  may  be  replaced  by 

z^l2x. 


The  corroction  on  E  ,  which  is  equivalent  to 
n 

«  a  ) 

Ti  n  n 

was  derived  fron  a  theoretical  consideration  of  the  velocity  profile  downstream 
of  a  wire.  This  was  rado  alon?  the  linos  presented  by  S.  I,  Pai  (Hef.  7)  for 
two-dinonsional  jot  nixing,  replacing  tlie  jet  in  a  uniforn  stroan  by  an  obstacle 
(vnre)  in  a  uniforn  stream.  Tiie  constant  A*  in  the  above  equation  is  a  function 
cf  Uio  Uiernal  conductivity  and  temperature  of  the  gas,  the  heat  release  rate  of 
the  wire  and  the  drag  on  tiie  wire.  Since  the  drag  on  the  wire  is  not  known.  A' 
(thus  a)  is  determined  by  requiring  A  to  be  tlie  sane  at  two  separated  x  values. 

References 

1)  Hirschfelder,  J.  0.,  Curtiss,  C.  F.,  and  Bird,  R.  B.,  The  Molecular  Tneory 
of  Gases  and  Liquids.  (John  Wiley  &  Sons,  Inc.,  New  York,  195M. 

2)  V/alker,  R.  E.  and  V/estenberg,  A.  A.,  J.  Chem.  Phys.  1139  (1958). 

3)  Task  R  Quarterly  Progress  Rejwrt  No.  2,  1  July  *■  30  September  1959,  The 
Johns  Hopkins  University,  Applied  Physics  Laboratory  Report  No.  TG-331-2. 

h)  Task  R  Quarterly  Progress  Reoort  No.  1,  1  April  -  30  June  1959,  The  Johns 
Hopkins  University,  Applied  Physics  Laboratory  Report  No.  TG-331-1. 

5)  National  Bureau  of  Standards  Circular  561i,  Tables  of  Thermal  Properties  of 
Gases, 

6)  Johnston,  ’i.  L.  and  Grilly,  E,  R.^J.  Chem.  Phys.  ]^,  233  (I9l6). 

7)  Pai,  S.  I,,  J.  Aero  Sci.  16,  16}  (1919), 


-  19  - 


in.  TRA!SPJKT  PHOr-i^TY  STUDIES  IH  DISSCCIAT-D  GASES 


(R,  E,  Vfalkor  ai>d  A.  A.  -'testcnbare) 


Ob:ioctivc 

ConvoctivQ  heal  transfer  from  dissociated  ^asos  to  a  catal^'tic  vail  is 
significantly  affected  by  the  clienical  kinetics,  the  nature  of  the  catalytic 
surface,  and  the  transport  properties  of  the  rdxture.  Although  gas  nJiase 
kinetics  for  dissociated  gases  have  boon  studied  separately  in  shock  tube 
exporinents,  attempts  to  neasurc  the  catalytic  nature  of  a  sui'face  liave  re¬ 
vealed  that  it  is  difficult  to  separate  the  catalytic  effects  from  those  due 
to  diffusion.  Until  very  recently  (Ref,  1)  no  separate  measurements  of  dif¬ 
fusion  coefficients  had  been  made,  so  that  gross  estimates  of  this  transport 
property  are  generally  required  to  obtain  catalytic  activity  of  various  surfaces. 
Additional  diffusion  measurements  for  a  variety  of  dissociated  gases  are  de¬ 
sirable  for  this  reason. 

In  addition,  measurements  of  diffusion  coefficients  as  a  function  of 
temperature  is  an  excellent  way  of  determining  intermolecular  forces  between 
unlike  molecules.  For  dissociated  gases,  the  (labile)  atom-(parent)  molecule 
intermolecular  forces  could  be  obtained  in  this  manner  and  subsequently  be  used 
to  predict  other  related  transport  coefficients,  (These  intermolecular  forces 
can  be  determined  from  first  principles  for  only  the  simplest  interactions 
(H-H,  H~He).)  Measurements  of  thermal  conductivity  and/or  heat  transfer 

coefficients  in  dissociated  gases  are  also  of  considerable  intrinsic  and  prac¬ 
tical  interest.  The  present  project  is  aimed  at  exploring  and  using  various 
methods  of  experimental  measurement  of  transport  properties  in  dissociated  gases. 
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Exnorinental  Todmiguc 

£x|»riTO;.t^  effort  thus  far  lias  been  restricted  to  elementary  experi¬ 
ments  with  dissociiiteh  gases  created  by  a  radio  frequency  (HI*)  discitarge.  A 
discussion  of  tliose  jirollnina.ry  exporinents  is  {jiven  below, 

Pi*olicdnary  axcerinonts 

In  order  to  obtain  a  "fetU."  for  tlie  concentration  of  atoms  t'nat  can  be 
produced  by  an  HF  discliarge  and  the  lifctinos  of  tliese  atoms,  the  sinple  flow 
system  shown  in  Fig.  III-l  lias  been  used.  Atoms  are  generated  by  passing  a 
diatonic  gas  (O2  and  II2  have  been  used)  through  the  inductance  of  a  tuned  tank 
circuit  coupled  to  an  RF  transmitter  capable  of  yO  watts  input  to  final  stage 
(EICO  #720).*  Atom  are  detected  downstream  of  the  discharge  by  a  fine  platinum 
(c.OOl  in.  dia.)  maintained  at  a  fixed  temperature  by  a  variable  healing 
current.  The  change  in  heat  dissipation  with  and  without  the  discharge  is  a 
measure  of  heat  release  through  catalysis  (atom  recombination  on  the  surface 
of  the  wire)  and  the  changes  in  convective  and  conductive  heat  transfer.  The 
convective  and  conductive  heat  transfer  is  eliminated  by  using  a  similarly 
heated  platinum  wire  that  is  coated  with  silica  to  render  it  non-catalytic. 

The  heat  release  through  catalysis,  Q,  can  Uien  be  related  to  the  number  ol‘ 
atoms  by  the  expression 

Q  =  ^nca)h  (1) 

8kT\^^^ 

where  n  is  the  number  of  atoms  per  unit  volume,  c  *=  ( - )  is  the  mean 

TT  m 

molecular  speed,  a  is  the  surface  area  of  the  wire,  y  is  the  surface  catalytic 
activity  for  atom  recombination  (i.e.,  the  percentage  of  atoms  striking  the 
surface  which  recombine)  and  h  is  the  dissociation  energy  per  molecule, 

A  total  of  500  watts  is  now  available  by  the  addition  of  a  recently  purchased 
RF  amplifier  (Viking  Courier), 
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Prelinianry  oxp9rin9ntn  have  been  restricted  to  oxygen  gas  and  the  effects 
of  pressure  at  different  flow  rates.  A  sua’.Tai7  of  Uw  data  is  given  in  Fig.  111-2. 
At  pressures  below  about  0,J.  m  !!g  the  efficiency  of  the  .’tF  disdiarge  is  quite 
poor.  At  the  elevated  pressures  the  discimrge  is  also  less  efficienti  and  in 
additioni  atora  are  presuraably  lost  by  wall  rccorsbination.  Additional  refined 
experiiaents  along  this  line  are  planned. 

Proposed  Exporlnent  for  Diffusion  Hoasurenents  of  Dissociated  Species 

The  usual  nothod  of  neasuring  the  efficiency  of  a  surface  for  aton  recon- 
bination  is  to  place  the  catalytic  surface  at  the  closed  end  of  a  side-arm 
attached  to  a  tube  containing  a  flowing  dissociated  gas  (itef.  2).  The  net  flow 
in  the  side  arm  is  zero  and  the  flux  of  labile  species  to  the  catalytic  surface 
is  controlled  by  the  diffusion  through  the  gas  mixture.  By  measuring  the  heat 
release  at  the  surface  (via  a  temperature  rise)  a  measure  of  the  flux  is  obtained, 
and  by  varying  the  depth  of  the  side  arm  the  ratio  of  catalytic  activity  to  dif¬ 
fusion  coefficient  can  be  determined.  A  possible  refinement  of  this  method  is 
to  alter  the  atom  flux  by  imposing  a  counterflow  (see  Fig.  III-3).  The  concen¬ 
tration  of  atoms  is  then  given  by 


C  + 


C  Yc 

0 

TT 


c  /c\ 
0  1 


-^(L-Z) 


(2) 


where  and  Cj^  are  atom  mole  fractions  at  z  =  0  and  z  =  L,  respectively,  U  is 
the  countercurrent  gas  velocity,  and  D  is  the  diffusion  coefficient.  In  deriving 
Eq.  (2)  loss  of  atoms  to  the  wall  of  the  tube  by  heterogeneous  reactions  and  by 
gas-phase  reaction  has  been  neglected  and  it  is  assumed  that  C-=^-=^l.  For  small 
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-3  SKETCH  OF  SYSTEM  USING  HETKHOGENEOUS  HKACTION  WITH  COUSTKIlFLOW  AS 
METHOD  OF  MEASURING  DIFFUSION  COEFFICIENTS  OF  DISSOCIATED  GASES. 
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values  of  UL/D  £q.  (2)  expands  to  give 

(3) 


(h) 

although  this  condition  need  not  necessarily  apply.  By  measuring  a  quantity  that 
is  proportional  to  the  concentration  of  atone  at  the  surface  (say^  a  temperature 
rise)  and  by  accurately  metering  the  countercurrent  flow,  the  diffusion  coefficient 
can  be  measured. 

At  the  low  flow  rates  encountered  in  an  experiment  of  this  type,  the  flow 
nay  consist  entirely  of  boundary  layer  and  the  simple  mathematics  surely  will  not 
apply.  To  avoid  boundary  layer  problems,  it  nay  be  desirable  to  use  spherical 
geometry  for  which  an  expression  similar  to  Eq.  (3)  can  be  developed.  Both  geo¬ 
metrical  arrangements  are  presently  receiving  attention  as  possible  methods  of 
measuring  diffusion  coefficients  of  dissociated  gases. 
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For  some  cases,  so  that  Eq.  (3)  reduces  to 
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Objectivo 

One  of  the  difficulties  in  the  way  oi  adequate  understanding  of  U)o  aerc- 
dynamics  of  high  tonperaturo,  high  speed  gas  flow  has  been  Uie  lack  of  inforpsation 
about  the  chemical  kinetics  involved.  Obtaining  this  inloroation  has  been  dif¬ 
ficult,  especially  in  those  cases  whore  premixino  the  reactants  before  bringing 
then  to  reaction  temperature  is  impossible  because  of  the  rapidity  with  which 
the  reaction  proceeds.  Laminar  flame  studies  such  as  are  discussed  in  Section  I, 
are  one  way  of  acquiring;  high  ttr,;>?rature  ki’^eiic  data  if  t:io  reactants  can  be 
premixed.  Another  approac;:  V  '  adapi  ari-  tia-  point  source  teonr.ique 

(Ref,  1)  used  for  the  measurement  of  ciiffusion  coefficients,  a  method  wliich  should 
be  especially  useful  when  premixing  is  impossible.  This  technique  will  allow 
preheating  of  the  reactants  before  mixing,  steady-state  operation,  a  fair  tem¬ 
perature  range,  and  simple  psuedo-first  order  kinetics  for  analysis.  Thus  it 
can  be  seen  that  information  obtained  by  this  method  might  be  a  useful  supplement 
to  that  obtained  in  the  laminar  flame  work. 

Theory  of  the  Method 

Figure  IV -1  is  a  schematic  representation  of  the  principle  it  is  proposed 
to  use.  Here  a  trace  amount  of  one  reactant  is  continuously  injected  into  the 
center  of  a  laminar  stream  of  a  second  reactant.  Under  conditions  of  constant 
temperature,  pressure,  and  uniform  carrier  gas  velocity  U,  and  when  the  mixture 
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SCHEMATIC  DIAGRAM  OF  REACTING  POINT  SOURCE  SYSTEM 
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is  so  dlluto  in  trace  gas  that  the  ticnsiiy  is  ossenliaUv  constant,  the 
equation  governing  the  cerxcntrution  C  (rsolo  lYnction)  of  trace  gas  is 


where  D  is  the  binary  diffusion  coefficient  (constant)  and  k  is  the  pseudo 
first  order  rate  constant  for  the  disappearance  of  trace  gas,  i,e.,  the 
chemical  rate  expression  is  taken  to  be 


dC 

3t  ' 


-  k  C 


(2) 


in  view  of  the  large  excess  of  carrier.  At  a  given  toaperaturc,  Eq,  (1)  ray 
be  solved  to  give 
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vjhere  Q  is  the  volumetric  flow  rate  of  trace  gas  from  the  source. 

To  obtain  k,  one  would  measure  C  first  at  a  series  of  position  on  the 
axis  where  x  =  r,  and  second,  at  a  series  of  positions  r  along  a  transverse 
profile  of  fixed  x. 

For  the  first  case,  one  may  shov/  that 


k  = 


il) 


vjhere  is  the  slope  of  the  line  obtained  from  plotting 
^2  ^2 

log -  versus  (xp  -  x^)  ^ 


where  C  =  for  x  =  x^ 
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For  tho  second  case, 

2  ^ 


C<>r, 


k  -  «  U-AD 

where  slope  obtained  by  plottii^  log  versus  (r  >r  ) 

1  2^  * 

where  C  ■  for  r  ■  r^^  . 


(5) 


Thus,  if  Q  and  U  are  unknown  but  consUnt,  it  is  possible  to  dotemino 
tho  rate  constant  k  if  tho  dilfusion  coefficient  has  been  previously  doteminod. 
An  important  point  is  Wiat  only  ratios  of  concentrations  are  required,  so  no 
calibrations  are  necessary. 

Experimental 

One  system  of  interest  in  combustion  is  the  oxidation  of  carbon  monoxide. 
Inasmuch  as  the  diffusion  coefficient  of  CO  into  oxygen  had  already  been 
determined  with  precision  (Ref.  2)  and  some  rate  data  were  available  from  laminar 
flame  studies,  this  pair  was  chosen  for  the  initial  work, 

A  thermal  conductivity  cell  is  used  to  analyze  the  gas  mixture.  However, 
the  thermal  conductivity  of  CC  is  too  close  to  that  of  0^  for  the  method  to  be 
very  accurate.  Consequently,  the  follov/ing  scheme  was  devised.  The  gas  mixture 
drawn  into  the  probe  is  first  passed  through  a  liquid  nitrogen  trap  to  remove 
any  carbon  dioxide  which  has  been  formed.  The  mixture  is  then  passed  through 
a  heated  platinum  catalyst  to  oxidize  the  previously  unreacted  CO,  and  this 
mixture  of  CO^  and  0^  goes  into  one  side  of  the  T-C  cell.  The  gas  then  passes 
through  a  second  liquid  nitrogen  trap,  removing  the  CO2  and  allowing  pure  oxygen, 
the  reference  gas,  to  flow  into  the  other  side  of  the  T-C  cell.  This  continuous 
flow  system  thus  gives  a  record  of  the  CO  concentration  at  any  given  time. 
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It  la  apparorti  that  the  reacUni;  3a:iplc  drawn  Into  the  probe  at  ai\y 
point  RUst  be  quenched  offoctlvoly  In  order  for  neanineful  Roaaurononls  to 
be  nade.  The  aanpllnc  probes  used  will  be  slnilar  to  these  used  in  flsHi 

studios  (Rof.  3)y  where  sonic  flow  and  rapid  drop  In  pressure  on  entering 

the  probe  are  relied  on  to  quench  reactions, 

Rosulta 

As  of  the  present  tlno>  the  equipnont  has  been  assorhled  and  ciost  of 
the  operating  difficulties  ellninated.  However,  a  feu  still  remain  and  so  no 
results  of  ai\y  consequence  have  yet  been  obtained. 
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V.  RCCKS7  KOZZLE  FLUID  DY^^AIHCS 
(r,  K,  Hilly  H,  A.  V/allskoCy  <ind  H.  J,  U;icer) 

Objoctlve 

The  invosti cations  carried  out  in  this  problem  are  concerned  with 
the  dynamics  of  the  liich  swed  cas  flows  occurrinc  in  solid  pi*opollar.t 
rocket  nozzles.  V/ith  ivigh  temperatures  and  laree  mass  flow  rates  and  with 
rapid  expansion  of  the  gas  through  the  nozzle,  it  is  Important  to  describe 
the  fluid  properties  taking  into  account  non-equilibrium  and  non-ideal  gas 
conditions.  Experimental  data  pertinent  to  an  accurate  description  of  the 
gas  are  the  first  objective  so  that  an  adequate  theory  of  the  fluid  dynamic 
behavior  may  be  formulated. 

In  order  to  pro’/ide  realistic  and  representative  testing  conditions 
for  the  experiment  a  double  -base  solid  proiajllant  (AltP)*  has  been  chosen 
which  provides  the  five  basic  coratittients  common  to  most  high  impulse 
propellants.  These  are  hytirogen,  nitrogen,  carbon  monoxide,  carbon 
dioxide  and  water  vapor.  Impurities  of  the  order  of  ?  to  Jh  are  present 
as  is  the  case  in  all  propel' ants  and  it  is  possible  to  add  in  well  con¬ 
trolled  amounts  and  sizes  solid  particles  for  future  extensions  of  the 


The  propellant  is  manufactured  and  proof  tested  at  Allegany  Ballistics 
Laboratory. 


V 
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experlnental  studios.  To  begin  with  U  is  beliovod  advisable  to  keep 
the  problenis  as  uncorapli eaten  as  possible  wSiile  still  retaining  the 
fundanental  aspects  of  the  p!w!no?»na  under  investigation,  Costbustion 
pressure  and  temperature  of  Ute  grain  are  ICOJ  psi,  noninal  and  2500°K, 
respectively;  end  burning  is  cnplo.'.'cd  fee*  cither  10-soconci  or  30'-nocond 
operation.  These  conditions  are  sufficient  to  introduce  measurable 
effects  due  to  the  variable  gas  properties  as  the  gas  expands  Uu-ougii 
the  nozzle,  gas  non-equilibriun  conditions  in  the  expanded  flow  and 
heat  transfer  and  nozzle  divergence  angle  effects.  Associated  plienomena, 
such  as  erosion  and  deposition  of  solid  oarticles,  are  also  present  for 
observations  providing  additiotul  data. 

It  is  anticipated  that  as  the  v;ork  progresses  higher  specific 
impulse  propellants  vri.ll  be  utilized  providing  temperatures  in  the 
3000°K  to  liOOO°K  range,  and  that  investigatiens  of  the  effects  of 
metal  additives  to  the  grain  vrill  be  carried  out.  Cooling  and  heat 
transfer  measurements  to  the  nozzle  walls  will  be  irade  and  some  pro¬ 
perties  of  materials  directly  associated  vrith  the  rocket  gases  may  be 


undertaken 
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Introduction 

As  described  In  earlier  reports  (Refs.  1,2),  ore  of  the  problera  to  be 
undertaken  vas  the  determination  of  U\q  actual  flow  properties  in  a  rocket 
noetle  and  their  comparison  with  the  calculated  properties  based  on  gas  equi¬ 
librium  assumptions.  It  vas  evident  Uiat  with  high  inpulsc  propellants  which, 
during  burning^ produce  a  high  enUmlpy  gas  nixture,  tiie  rate  of  conversion 
of  this  thermal  energy  to  kinetic  onerg>'  must  govern  the  over-all  efficiency 
as  well  as  the  environmental  conditions  of  the  nozzle  and  associated  equipment. 
Any  non-equilibrium  condition,  e.g.,  the  relaxation  of  excited  states,  or  ai\}' 
chemical  reactions  which  occur  during  the  expansion  process  in  the  nozzle  nay 
be  significant  in  the  ultimate  kinetic  energy  of  the  exhaust  gases.  Inter¬ 
related  with  these  processes  are  surface  plienomcna  of  the  containing  walls, 
including  heat  transfer  and  possible  efficient  methods  of  cooling  and/or 
insulation. 

Nozzle  Surveys 

The  experimental  work  was  directed  toward  making  accurate  aerodynamic 
surveys  of  the  flow  Mach  number  and  pressure  distribution  in  the  nozzle  where 
both  the  gas  density  and  static  temperature  were  relatively  high.  Two  series 
of  tests  were  made  using  nozzles  of  12.5°  and  26°  total  divergence  angles,  so 
that  the  flow  or  transit  time  of  the  gas  from  throat  to  survey  station  was  varied 
by  a  factor  of  two,  whereas  Mach  number  and  flow  velocity  would  be  approximately 
the  same.  Because  of  the  severe  environmental  conditions,  fixed  probes  were  used 
and  not  all  of  the  data  obtained  were  valid  due  to  mechanical  failures.  However, 
although  one  or  two  tests  remain  for  the  12.5°  nozzle,  the  survey  in  both  nozzles 
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at  the  M  ■  3  station  are  sufficiently'  conplete  to  provido  a  significant  conparison 
The  data  are  shown  in  the  top  part  of  Fig.  V>1  with  one  set  for  Uio  12.5^  noccle 
and  two  sets  for  the  2(P  noszle.  The  tests  with  the  26°  nocslo  were  carried  out 
firsts  and  only  one  noszle  was  used  in  which  Uiero  was  a  progressive  enlargcnent 
of  the  throat  diameter  (  0.001"  per  10  second  run).  Co!isequently|  there  was  a 
slight  decrease  in  area  ratio  as  the  tests  progressed;  also,  from  ex;«rience 
Inprovemnts  were  nade  in  the  rccchanical  support  of  Uie  probes.  Tltc  line 
through  U)e  lower  data  represents  the  first  tests  and  tiio  upper  line  the  final 
tests,  which  are  believed  to  be  a  more  accurate  description  of  the  26°  nozzle 
flow.  Nevertheless,  because  a  comparison  between  the  two  nozzles  of  differing 
divergence  angle  is  to  be  made,  the  conservative  view  is  to  take  the  lower  curve 
for  the  26°  nozzle.  The  data  displayed  in  graphical  form  are  the  pitot- to-supply* 
pressure  ratios,  ^  function  of  distance  from  the  nozzle  axis.  Tlie  data 

shown  in  Table  I  are  obtained  from  the  wall  static  pressure  measurements  P^. 

TABLE  I 


Nozzle 

Test 

P  /P 

S  0 

H 

26° 

58 

0.0225 

2.99 

66 

0.021. 

2.95 

68 

0.023 

2.98 

12.5° 

76 

0.021 

2.95 

77 

0.025 

2.93 

78 

0.02ii 

2.95 

This  measurement  provides  the  expansion  ratio  and  enables  the  local  Mach  number 
to  be  determined.  The  Mach  numbers  as  tabulated  and  as  indicated  in  Fig.  V-1 
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aro  calculated  fron  tiio  prcosuro  ratios  assuni%'  tiiat  thQi*nody;ianic  (Hiuilibriittj 
obtains  throughout  Urn  expansion  process,  and  Umt  Uje  conposition  is  iVonon 
at  proportions  conputed  for  1500®K  and  is  in  local  equilibriun  at  higher  tcnpora- 
lures.  Witliin  the  accuracy  of  naasurenenl  the  Hach  nu.'^bers  and  expansion  ratios 
in  the  tuo  nozzles  are  comparable,  but  with  a  slight  trend  towaixi  higher  Macii 
numbers  for  tiie  larger  divergence  angle. 

In  order  to  compare  the  pitot  ;u:ossure  distribution  data,  the  data  are 
shown  in  the  bottom  part  of  Fig,  V-1  as  the  change  in  pitot  pressure  ratios  re¬ 
sulting  from  the  change  in  divergence  angle.  It  is  noted  that  diere  is  a  ciiange 
in  measured  pitot  pressure  of  nearly  6^  in  the  central  region  of  the  flow  for 
these  surveys  conducted  at  eo.ual  geometric  expansion  ratios.  The  lower  curve 
labeled  "ideal  conical  flow"  is  the  pitot  pressure  difference  that  would  be 
expected  on  the  basis  of  pure,  simple  conical  flow.  The  difference  between 
these  two  curves,  approximately  1,8^  centerline,  nay  be  attributed 

to  non-equilibrium  effects.  Had  the  upper  line  of  data  for  the  26°  nozzle  been 
used,  they  would  have  indicated  a  difference  over  the  conical  flow  effects  of 
approximately  %  at  the  centerline. 

It  may  be  noted  that  in  the  calculation  of  the  flow  properties  of  the 
gas  mixture  as  it  expands  through  the  nozzle,  the  composition  of  the  gas  has 
been  assumed  to  be  in  equilibrium  down  to  1500°K  and  frozen  thereafter.  The 
contribution  of  heat  released  (or  absorbed)  by  chemical  reactions  is  not  taken 
into  account.  In  this  particular  case,  if  the  gas  does  in  fact  maintain  equi¬ 
librium,  the  water  gas  equation 

HgO  +  CO  ^  +  CO2 


(1) 
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Fig.  V-1  TOP:  PRESSURE  SURVEYS  IN  ROCKET  NOZZLES  OF  26°  AND 

12.5°  DIVERGENCE  ANGLES. 

(EXPANSION  RATIO  =  5.85,  P^  =  1150  PSI,  T^  =  2500°K) . 

BOTTOM;  PITOT  PRESSURE  DIFFERENCE 
BETWEEN  26°  AND  12-5°  NOZZLES. 
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indicates  that  fron  tSie  throat  (tonporature-^  2250®K)  to  downstream  position 
wiscre  Uie  tonparaturc  is  15(X)°K,  a  heat  release  ol*  some  71  DTU/sec  would  occur, 
Ti)is  is  of  the  oi'dor  of  h.S%  of  t!®  ki:Totic  enorjiv  of  tJjo  strean,  V/Jietiier  or 
not  this  occurs  depends  on  the  reaction  rates.  Just  downs trean  of  tl)e  throat 
the  cas  velocity  is  approximately*  0.2  Inch/p  sec  and  at  the  su2*vey  station 
about  0.1/tsec.  Thus,  Uie  significant  tire  scale  is  of  the  order  of  10  sec, 
since  the  distance  travelled  by  the  gas  is  several  indies  from  tiiroat  to 
survey  station.  For  the  conditions  which  are  inposed  by  the  propellant  in 
this  case,  it  is  expected,  therefore,  that  effects  due  to  chemical  non- 
equilibrium  will  be  restricted  to  orders  of  a  few  per  cent.  With  higher 
energy  propellants  where  the  oonbustion  temperatures  range  up  to  bOOO^Kj 
it  is  expected  that  non-equilibrium  effects  nay  be  much  more  pronounced 
and  that  their  acciu-ale  determination  nay  be  of  importance  in  the  nozzle 
design,  in  the  materials  used  to  contain  the  gases,  and  in  the  heat  transfer 
properties. 

The  pressure  surveys  to  date  have  been  limited  to  a  single  transverse 
plane,  and  for  this  reason,  it  is  difficult  to  say  whether  the  measured  net 
difference  is  entirely  due  to  non-equilibrium  chemical  or  thermal  effects  or 
perhaps  due  to  different  degrees  of  distortion  from  ideal  conical  flow,  Ad- 
ditioaal  surveys  are  planned  at  two  more  stations  in  the  nozzle  as  well  as  an 
effort  tc  determine  any  flow  distortion  effects  of  throat  and  nozzle  contour 


in  the  subsonic  area. 


Gas  Composition 


Very  closely  connected  ifith  non-equilibrium  condUions  is  the  gas  conpo  • 
silion  in  toms  of  t‘10  concentration  of  tlw  various  species  present  and  tJieir 
state.  A  dotornination  of  the  composition  from  point  to  point  along  the  nozzle 
axis  requires  essentially  an  instantaneous  noasuroment,  so  t5je  convenlioiml 
canpling  tediniques  are  not  practical,  Spectrographic  techniques,  on  the  other 
?uind,  appear  pronising,  and  several  of  the  species  are  known  to  have  absorption 
Ixands  in  the  infrared.  As  a  result,  it  was  decided  to  use  a  scanning  sjx»c- 
trorjotor  to  attempt  a  quantitative  dotemination  of  the  concentration  of  those 
sjiecies  witii  absorption  bands  in  the  Ip  to  l5p  region.  The  instrumentation 
problem  was  believed  to  lie  mainly  in  accomodating  the  required  optical  and 
electric  equipment  to  the  enviroamental  conditions  and  in  developing  suitable 
techniques  for  quantitative  measurements.  A  high  airborne  sound  level,  large 
non-uniform  temperature  gradjents,  short  time  e>q)osure,  and  remote  operation 
constituted  the  environmental  nroblems.  Since  no  suitable  instrument  was 
available  commercially,  the  vfork  required  the  adaptation  and  modification  of 
basic  components  to  an  apparatus  v;hich  could  survive  the  conditions  imposed 
and  display  the  information  gathered  in  a  readily  usable  form.  Since  the  rocket 
gas  is  a  mLxture  of  M2,  H2>  CO,  CO25  and  H2O,  three  of  these  gases  should  pro¬ 
vide  measurable  absorption  bands  in  the  infrared.  Furthermore,  since  the 
relative  abundance  of  K2  does  not  change  during  the  expansion  process,  the 
quantitative  measurement  of  the  three  I-R  absorbers  would  lead  to  a  determi¬ 
nation  of  the  relative  amounts  of  all  five  gas  species. 

In  all,  nine  environmental  tests  have  been  run  with  the  scanning  spec¬ 
trometer  (Ref.  3)  built  at,  APL  several  years  ago.  There  have  been  a  succession 


of  nodifications  nado  as  the  tests  pro^jossed,  until  at  the  present  lino  iho 
uiTangement  of  equipaont  is  as  shown  scheraatically  in  rig.  V-2,  Tlw  earlier 
tests  did  not  obtain  spectra  of  the  rocket  gases,  bo  cause  they  wore  dosigtwd 
only  to  operate  the  equirrent  lookir^  at  a  duiro'  source  (fJernst  glower)  wfion 
tlie  systen  was  exposed  to  the  150  db  noise  level.  ?h.o  last  four  tests  provided 
3j)ectra  of  the  gas,  and  sanpleu  of  the  results  obtained  are  ,.ivon  in  Kig,  V-3, 

As  indicated  in  rig.  V-2,  a  ono-iulf  inch,  herds  pile  rically  tipped  cylinder  wodel 
was  heated  by  the  strean  to  incandescence,  and  then  erdtted  radiation  which 
traversed  the  pas  flovdng  around  it  and  mssed  tlirough  the  tunnel  window  to 
the  spectrometer.  Since  H2O  and  CO2  are  present  in  the  atrosphere,  Uie  optical 
path  for  test  section  to  spectrorneter  is  covered  with  a  tube  and  flushed  with 
dry  N2.  Similarly,  the  optical  system  of  the  spectrometer  is  sealed  and  flushed 
with  N^.  The  dual  radiation  pyrometer  also  monitored  the  radiation  from  the 
model  on  the  opposite  side,  so  that  a  time  history  of  its  temperature  was  ob¬ 
tained  and  is  shown  at  the  bottom  of  Fig.  V~3.  The  model  rises  in  temperature 
with  time  as  shown  in  the  temperature  record  and  by  the  magnitude  of  the  conti¬ 
nuous  spectral  background  as  shown  in  the  traces  for  boWi  T/N  71  and  72, 

The  upper  traces  in  Fig.  V-3  were  copied  from  polaroid  photographs  of 
the  C.R.O.  trace  during  Test  71,  and  while  these  traces  have  lost  in  definition 
due  to  the  copying  process,  the  originals  were  sufficiently  detailed  to  permit 
the  identification  of  the  bands  shown.  For  Test  72,  the  resolving  power  of  the 
spectrometer  was  increased  by  cutting  the  entrant  slit  width  in  half,  and  the 
results  are  shown  in  the  sketches  taken  from  35mm  movie  records  made  throughout 
the  run.  The  double  trace  seen  in  this  test  was  due  to  the  fact  that  the  camera 


SCHEMATIC  LAYOUT  OF  RAPID  SCAN  SPECTROMETER  FOR 
VIEWING  ROCKET  GASES 
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Fig.  V-3  SPECTROMETER  TRACES  RECORDED  FOR  TESTS  71  AND  12. 
TEMPERATURE  RECORD  OF  SOURCE  FOR  TEST  72. 


shutter  vas  operatij^  at  1?  frares/sec,  vimntis  scannlj^  rate  vss  6J  C}j@, 
Ttere  was  also  sone  ovldcjjce  oi*  rsicropSionics  in  tnees  traces,  which  nusi  I*© 
olinlmteil  before  Uw  absorption  bajsiis  can  be  clearly  resol veu.  It  sSwuld  be 
ronenbered  Uiat  the  total  rjnninj  tine  for  liwsQ  tests  was  10  seconds*  Tins 
quarts  tost  section  window  and  a  quarts  coUinatin^  lens  Unit  the  wave  len^U) 
to  about  3*5  M*  It  is  planned  to  replace  the  test  section  window  with  a  noro 
transparent  nediuA  in  the  infrared  (iCaCl  or  sapphire)  and  to  use  a  nirror  sy*sten 
for  collimtion*  In  order  to  provide  a  constant  enoro*  background  with  wave 
lengthy  it  is  also  planned  to  progran  the  slit  width.  Tiw  results  of  these 
tests  so  far  have 

a)  established  the  feasibility  of  this  absorption  technique, 

b)  indicated  several  inprovcnents  which  must  be  rado  for 
quantitative  data  on  the  oocposition,  and 

c)  shown  the  necessity  for  acoustic  shielding  to  reduce  micro- 
phonics  in  the  electronic  system,  whereas  toechanlcal  stability 
of  the  optical  system  seems  to  be  satisfactory. 

Fabrication 

During  this  period  the  new  rocket  combustion  chamber,  test  stands,  and 
adjustable  optical  benches  were  completed  and  installed  in  the  test  cell.  This 
equipment  v^as  designed  to  provide  a  greater  degree  of  safety  and  to  permit  the 
accurate  alignment  of  optical  equipment  for  the  tests.  These  new  combustion 
chambers  are  undergoing  proof  testing  at  Allegany  Ballistics  Lab  and  have  not 
been  used  as  .vet.  It  is  expected  that  the  new  equipment  will  considerably 
facilitate  the  preparation  and  conduct  of  the  testing  work. 
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NOTICE:  WHEN  GOVERNMENT  OR  OTHER  DRAWINGS,  SPECIFICATIONS  OR 
OTHER  DATA  ARE  USED  FOR  ANY  PURPOSE  OlHER  THAN  IN  CONNECTION 
WITH  A  DEFINITELrY  RELATED  GOVERNMENT  PROCUREMENT  OPERATION, 
THE  U.  S.  GOVERNMENT  THEREBY  INCURS  NO  RESPONSIBILITY,  NOR  ANY 
OBLIGATION  WHATSOEVER;  AND  THE  FACT  THAT  THE  GOVERNMENT  MAY 
HAVE  FORMULATED,  FURNISHED,  OR  IN  ANY  WAY  SUPPLIED  THE  SAID 
DRAWINGS,  SPECIFICATIONS,  OR  OTHER  DATA  IS  NOT  TO  BE  REGARDED  BY 
IMPLICATION  OR  OTHERWISE  AS  IN  ANY  MANNiSIR  LICENSING  THE  HOLDER 
OR  ANY  OTHER  PERSON  OR  CORPORATION,  OR  CONVEYING  ANY  ftiGHTS  OR 
PERMISSION  TO  MANUFACTURE,  USE  OR  SELL  AilY  PATENTED  INVENTION 
THAT  MAY  IN  ANY  WAY  BE  RELATED  THERETO.  * 


